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Abstract. In this paper we compute the mean field phase diagram of a doped 
antiferromagnet, in a magnetic field and with anisotropic exchange. We show that 
at zero temperature there is a metamagnetic transition from the antiferromagnetic 
configuration along the z direction to a spin-flop state. In the spin flop phase the 
system prefers a commensurate magnetic order, at low doping, whereas at larger 
doping the incommensurate phase is favorable. Contrary to the pure Heisenberg 
case, the spin flop region does not span an infinite area in the (A, h) plane, where 
A is the exchange anisotropy and h is the external magnetic field. We characterize 
the magnetic and charge-transport properties of the spin-flop phase, computing the 
magnetic susceptibility and the Drude weight. This latter quantity presents a sudden 
variation as the spin-flop to paramagnet phase transition line is crossed. This effect 
could be used as a possible source of large magneto-resistance. Our findings may have 
some relevance for doped La2-aSr,5Cu04 in a magnetic field. 
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1. Introduction 

Metamagnetic transitions are ubiquitous in Nature. We can find them in pure magnets, 
in boson systems, in spin density wave systems, and in doped antiferromagnets. 

When placed in a magnetic field many magnetic materials undergo first-order phase 
transitions. These materials are called metamagnets. Two of the best studied magnetic 
materials exhibiting metamagnetic transitions are MnF2 and FeCl2. El 0] In FeCl2 
there is a first-order transition from an antiferromagnetic (AF) to a paramagnetic (P) 
state, and in MnF 2 there is a first-order transition from an AF to a spin-flop (SF) state 
(see figure ^ panel a, for a schematic idea on the spin configurations). In addition to the 
metamagnetic transitions above mentioned, these materials also undergo second-order 
phase transitions. The way the second- and first-order transition lines touch each other 
introduces different types of critical points. In MnF 2 and FeCl2 there is a bicritical and 
a tri critical point, respectively. Some materials, as GdAlC>3, jl] can present different 
types of critical points depending on the orientation of the magnetic field. 
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From the theoretical point of view, and as early as 1936, Neel predicted a first 
order magnetic transition for an anisotropic antiferromagnet in a magnetic field. [3] 
This metamagnetic transition was named spin-flop transition. Improvements of Neel's 
results were achieved using Green's functions and Holstein-Primakov bosons. U\ 
Later, a scaling theory for bicritical points, based on a renormalization group analysis 
was introduced. [S] More recently, a zero-temperature spin-flop transition in square and 
cubic lattices was studied, using exact diagonalization, 9 and the SSE Monte Carlo 
method. [10J The finite temperature study of the model phase diagram, using exact 
(numerical) methods was done recently, but for the two dimensional case only, and in 
the context of hard core bosons. fT] (The three dimensional case, at finite temperatures, 
was also studied. [T^]) 

The exact mapping between hard-core bosons with nearest neighbor interaction 
and the spin 1/2 anisotropic Heisenberg model, JH] permits to obtain results for both 
physical systems, from the study of one of them alone. In the boson language, the 
antiferromagnetic and spin flop phases correspond to a Mott insulator phase (a solid 
phase), where the bosons are locked at the lattice sites by the nearest neighbor repulsion, 
and to the super-fluid phase, respectively. The observation of the Mott insulator to 
super-fluid transition (and vice versa) in a Bose-Einstein condensate of 87 Rb atoms 
confined by an optical trap,jT3] gave an enormous boost to the research, both theoretical 
and experimental, in this field. ^B] The effect of disorder in the phase diagram of the 
two dimensional boson-Hubbard system has been studied as well. [17 

Also in doped magnetic materials, such as in L^^Sr^CuO^, metamagnetic 
transitions have been observed. [T9"] Both in pure La 2 Cu0 4 and in lightly doped 
L^-aSr^CuO^ the spin flop transition has its origin in the Dzyaloshinskii-Moriya 
interaction. Although some theoretical study of this transition has been done in the past 
for pure La 2 Cu0 4 ,|2nj its study in doped L^-^Sr^CuO^ has not been pursued to our 
knowledge. The spin wave spectrum for doped L^CuO^ was consider Ivanov et oi..|21j 
This motivated us to carry general studies of metamagnetic transitions in anisotropic 
doped antiferromagnets. To the best of our knowledge, the only work dealing with 
metamagnetic transitions in strong-correlated electrons was done in the context of the 
Hubbard model. |22j 

The paper is organized as follows: in section El we introduce the model and the 
possible spin phases that may exist in the system; in section the zero- and finite- 
temperature phase- diagrams of the system are introduced and discussed, together with 
the effect of doping on the critical fields and on the spin-flop angle; in section 0] 
the metallic SF phase is characterized both in terms of their magnetic and transport 
properties. 
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2. Model and mean field equations 

Our starting point is a simple generalization of the usual t — t' — J model in three 

dimensions, where an anisotropic exchange term, of strength J A, in included 

t t' 
H = -- ( C *> c *+/V + H.c.) - - (Ci,*Ci+i3',v + H.c.) 

i,f3,tj i,f3' ,a 

+^ Y,( as i s Up + s?s? + p + sysip) - /i E 5? - a* E c t<v • (!) 

i,f3 i i,cr 

In the Hamiltonian (0), (3 and f3' stand for the lattice vectors connecting all nearest and 
second-nearest neighbors sites of site i in a 3D cubic lattice, respectively; the parameters 
t and t' stand for the hopping integrals connecting all nearest and second-nearest 
neighbors sites, respectively, h stands for the magnetic field and /i for the chemical 
potential. The density of electrons is defined as n = N e /N, where N e is the total 
number of electrons and N is the number of lattice sites. Since our calculations will 
be performed at electronic densities close to n — 1 it is convenient for latter use to 
introduce the doping of holes as 5=1— n < 1. 

It is a characteristic of t — J models (we are referring here to class of models) that 
only one electron at the most can exist at each lattice site. We enforce this constraint 
using a slave boson representation for the q j(J operators. [2S] In this representation we 
have c iiCT = b\f ij(J and S = \ J2a',a fi a ^a',afi,a', where a = (a x , a y , a z ) are the usual Pauli 
matrices. At each lattice site the total number of bosons and fermions is equal to one, 
and we have the constraint J2a fl afi,<j + = 1- In what follows, we shall consider that 
all bosons have condensed into the lowest energy state, and introduce (b) = Vo, where 
5 corresponds, for zero and moderate temperatures (see finite temperature discussion), 
to the doping introduced above. 

In order to introduce a mean field Hamiltonian suitable for the study of 
metamagnetic transitions we consider two sub-lattices, A and B, such that in each 
sub-lattice the spins are oriented as shown in Fig. ^ panel c. Referring to Fig. Q 
panel c, 9 = <fi = describes the AF configuration; 9 = n/2 and 7^ describes the 
SF configuration; 9 = 7r/2 and 4> = 7r/2 describes the P configuration (the term "P 
configuration" to name a state where both spins point up may cause some confusion, 
but in this context it means the system does not present a ferromagnetic ground state, 
this is, it does not possess a ferromagnetic order parameter, in the sense of Landau 
theory of phase transitions; the magnetic moment presented by the system is induced 
by the magnetic field); 9 7^ and 7^ corresponds to a mixed phase, where there is 
staggered magnetization in both the z and y directions. The spin averages of all these 
configurations, in sub-lattices A and B, are given by 

(S iA ) = (0, -S A sin(# -<f>), S A cos(9 - 0)) , 

iB ) = (0, S B sin(# + 0), -S B cos(# + 0)). (2) 

Note that we consider the magnitude of the average value of the spin to be different at 
the two sub lattices. This is needed because the magnetic field induces a certain amount 
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of ferrimagnetism in the system. The above phases are commensurate with the lattice, 
but it is well known that the isotropic t — J model supports spiral order in the xy plane 
with a momentum vector Q incommensurate with the lattice. ^H] To account for this 
possibility in the presence of magnetic field, we introduce averages of the spin operators, 
presenting spiral order in the xy plane 

(Si) = S (sm (j) cos 9i, sin sin 6^, cos0) , (3) 

where Q± = Q ■ Ri, and represents the angle of S with the z axis (see figure ^ panel 
b). This phase competes with the commensurate SF order. It is also possible that the 
commensurate AF order in the z directions may compete with an AF order presenting 
incommensurate spiral order in the xy plane 

(S iA ) = Sa (sin 4> A cos 9i, sin <p A sin 9 h cos 0a) , (4) 
(S il3 ) = S B (sin <p B cos 8i, sin <p B sin 6i, -cos0b). (5) 

In our study we have not found solutions for the mixed state and for the AF order 
presenting incommensurate spiral order in the xy plane. 

2.1. Commensurate AF and SF phases 

Using the averages (J2J) and introducing an Hartree-Fock decoupling of the Hamiltonian 
(JTJ we obtain, after Fourier transforming the operators, the following mean field 
Hamiltonian 



H M F = E e lO)(4,A<x + b\,a a k,a) + X)[ £ 2( fc ) + ah B\A,c a k,a + 
k,a k,cr 

5> 2 (fc) + ^M,a* + E frfl4, ff ak.-«T + E kA,*K-° ■ ( 6 ) 

k,a k,<r k,<r 

where the a\ a and the b ka operators refer to the sub-lattices A and B, respectively, k 
stands for k, the k summation runs over the magnetic Brillouin zone, and 

£i(k) = — 2tS cos ki , €2 (k) — —At' 8 ^ cos ki cos kj , 

i=x,y,z j^i=x,y,z 

h Z A/B — 3JA(Sa/ b ) — — , h x A / B = 3J (S A / B ) . 

We point out that in €\(k) and 62(h) the condensation of the bosons renormalizes the 
hoping integrals to t5 and t'5. Since we assumed the condensation of the bosons, the 
constraint (J2 a fiafi,a + b\bi = 1) gives an equation for the number of particles in terms 
of the doping S 

^((«U,) + (tw) = i-^ (7) 

k,cr 

In the P phase there is one mean field parameter only: the average value of the 
spin S = S A = S B ; in the AF and SF phases we have the mean field parameters S A and 
S B , and S and 0, respectively. All physical quantities characterizing model © can be 
obtained from the associated single particle Green's functions. 
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Figure 1. Panel a: schematic arrangement of the spin on the two sub lattices in the 
AF, mixed, and SF states. Panel b: definition of the angles for the spins in the spin 
flop incommensurate spiral state given by Eq. 10. Panel c: definitions of the angles 
for the spins, including the AF, the SF, and the mixed states. The letters Sa and Sb 
denote the spins on sub-lattices A and B respectively. The spin configuration refers to 
the commensurate states given by Eq. (j2J. 



2.2. Incommensurate (spiral) SF phase 

Considering the spin-flop incommensurate states given by 0, the mean field 
Hamiltonian can be cast in the form 



where 



H = 5>(fc + Q)fl +QA f k+QA + e^flj^] 

k 

+ A ^2(fk+Qjfk,l + f k ,ifk+Q,d + Eo , 



(8) 



e T (/c + Q) 


= ei(k + Q) + e 2 (k + Q)+h z -^, 


h z 


h J a n 

= - - + - AS cos (f>z, 




= ex(k) +e 2 (k) - h z - n , 


A 


= ^Ssin07(Q), 


l{k) 


= 2 X! cos(fcj) , 



i=x,y,z 
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E = -^AS 2 cos 2 (pzN-^S 2 sin 2 (P-f(Q)N. 

The mean field parameters are the amplitude S, the angle <fi and the incommensurate 
momentum Q = (Q x ,Qy,Qz)] z is the coordination number. The corresponding saddle 
point equations are obtained from the free energy, determined from 

^=-TEE Ml + e- E « {k)/T ) + Ml - S)N + E , (9) 

k a=± 

with E a equal to 

B. = M * + Q 2 )+e ' W + |,/h(* + «)-« l (*)P + 4i' (10) 

and are given in the appendix (the incommensurate momentum Q = (Q x , Q y , Q z ) is 
also determined from a saddle point equation). 

We note that incommensurate spiral order in the xy plane leads to a well defined 
mean field theory, since it's possible to define a close set of equations of motion for the 
Green's functions. Conversely, this is not possible for incommensurate states in the zx 
or zy planes (in a magnetic field). Due to symmetry, the incommensurate states can 
only be of the form (Q,Q,Q), (Q,Q,tt), and (Q, 7r, 7r). In Table (JTJ) we show the effect 
of the doping 5 on the incommensurate (Q, Q, Q) and (Q, n, n) wave vectors. It is clear 
from table [T] that the free energy for the three states (Q, Q, Q), (Q, Q, 7r), and (Q, n, it) 
should be essentially the same, and therefore the value of ^(Q, Q,n) is not presented. 
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Figure 2. Effect of doping and A on the value of the incommensurate momentum, for 
T = 0.002 (temperature in units of t). The t, t' , and J parameters are those of Table 
[TJ and in the left panel h = 0.2, 0.6 and A = 2; in the right panel n = 0.92 (5 = 0.08) 
and h—0.2, 0.3, 0.6. Please note that the scanning on S for h = 0.2, corresponds to 
the second column of tabled 
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Table 1. Free energy values for incommensurate momentum Q as function of density, 
for T — 0.002, h = 0.2, A = 2. Here and in the remaining figures we have used t = 1, 
t' = 0.1 and J = 0.1. These results have been obtained in a lattice of 100x100 xlOO. 



It is clear from table ^that the incommensurate spiral phase has a lower free energy 
for moderate doping. As the doping is reduced the system finds the commensurate phase 
energetically favorable. [23 ESJ A detailed evolution, as function of doping and A, of 
the incommensurate spiral phase for the momentum (Q, it, tt) is presented in figure|21 In 
this figure we present results for different values of the magnetic field (h =0.2, 0.3, 0.6). 
The behavior of Q with A is seen to be non monotonous. For a given value of h, tt — Q 
may present a minimum (Q maximum) for a given value of A. If h is large (for example, 
h = 0.6), such that as A is reduced the SP to P transition-line is crossed, we see Q 
tends to flow away from the commensurate tt value. In the other cases (for example, 
h = 0.3), Q approaches the commensurate value first (tt — Q is reduced) , reaches the 
value where tt — Q is minimum, and starts to deviate again from the commensurate 
value (tt — Q increases) . (In the P phase there is no meaning for the incommensurate 
or commensurate states.) 



3. Phase diagram at zero and finite temperatures 

At T = and 5 = it is a simple exercise to obtain an analytical solution for the 
phase diagram of the system. There is no spiral order and the paramagnetic (Ep), 
antiferromagnetic (Eaf), and spin flop (E$f) energies are given (for spin 1/2) by 

E P = ^JzA-h, Eaf — > (U) 

and 

Jz h 2 

EsF = "T " J,(l + A) ' (12) 

respectively. At zero temperature the average value of S, Sa, and Sp is 1/2. From results 
fTTj) and fP2*|) the zero temperature phase diagram of the system can be obtained. The 
line separating the AF and SF phases is given by h — zJ^A 2 — 1/2, and the SF and P 
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Figure 3. Effect of doping on the magnetic structure at zero temperature. The t, t', 
and J parameters are those of Table and A = 1.2 and h=0.l. For these values of 
the parameters the AF phase is the more stable phase. 

phases are separated by the line h = zJ(A + l)/2. When we dope the antiferromagnet, 
there is a reduction of the average values of S, Sa, and Sb relative to 1/2, even at 
zero temperature, as can be seen from figure El This behavior is well known for the 2D 
isotropic t — t' — J model, and has been used to explain, at the mean field level, the 
reduction of the Neel temperature in high-temperature superconductors, upon doping. 
Here, the presence of the magnetic field combined with doping forces the AF phase to 
a ferrimagnetic phase, where the average value of spin in the two sub-lattices is not the 
same even at zero temperature. This is in contrast to the pure Heisenberg case, where 
Sa = Sb = 1/2 at zero temperature in the presence of h. Upon increasing doping, the 
average values of S, Sa, and Sb are reduced to the paramagnetic value, fixed by the 
magnetic field. Above a given value of 5, there is no magnetic order in the ground state 
of the system, and its behavior it that of a collection of fully polarized independent 
electrons. 

This reduction of S, Sa, and Sb together with the combined effect of A and h has 
consequences for the zero temperature phase diagram of the doped antiferromagnet, 
when compared with the pure Heisenberg case discussed above. The picture is the 
following: doping introduces holes here and there in the lattice and as consequence 
some of the magnetic interactions due to the Heisenberg term cannot be fulfilled; as a 
consequence, and for a given A > 1, as the magnetic field increases, the system finds 
favorable to take advantage of the Zeeman energy and therefore the first order AF-SF 
transition should occur at a lower value of h; for the same reason, the SF phase cannot 
fulfill all the antiferromagnetic interactions in the xy plane, and therefore the second 
order SF-P transition should also occur at lower values of h. 

This picture is confirmed by the phase diagram of figure 0] In this figure, and for 
A roughly in the range 1 < A < 4 we see that both the first order AF-SF and the 
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Figure 4. Left panel: zero temperature phase diagram of the 3D t — t' — J model, 
for 1 — S = n = 0.92. The t, t' , and J parameters are those of Tabled The black 
box marks the point where the first-order AF to SF and the second-order SF to PM 
transition lines meet. For values of A above the black box there is a AF to PM second- 
order transition line only. The pure Hciscnbcrg phase diagram is represented by the 
dashed lines. The full lines represent the phase diagram for the doped case. Right 
panel: variation of the SF angle cf>, in units of 7r, along the AF to SF transition-line. 
The inset shows the phase diagram close to A = 1. A minute reentrant behaviour is 
seen. 

second order SF-P transition lines are pulled down to lower values of h relative to the 
pure Heisenberg case. Furthermore, and in contrast to the pure Heisenberg case, the 
SF phase does not span an infinite area in the (A, h) plane. There is a point where the 
second order SP-P line meets the first order AF-SF line (represented by a black square 
in figure EJ). From this point on we are left with an AF-P second order transition; below 
this line the spins in the AF phase are fully aligned with the field, but have different 
magnitudes for the two sub-lattices. For values of A < 1 the system finds preferable to 
be in de SF phase. We note that in figure 0] the spin flop phase is of incommensurate 
type, in agreement with the results of figure El In fact using figure El in connection with 
the phase diagram of figure 0] it is possible to see the evolution of the incommensurate 
momentum Q over the phase diagram. 

We can now ask, how does the transition lines change with the doping? This 
question has experimental relevance in connection with measured spin flop fields in 
La2- x Sr x Cu04,[in] where it was found that the critical field for the SF transition is 
reduced upon doping. In figure |S] we show both the effect of 5 on the SF-P line, for 
small values of A (left panel), and on the AF — SF for two values of A. We see 
that for values of the exchange A = 2 there is a reduction of h c upon doping (we 
remark this behavior is not restricted to the single value A = 2, but it exists over a 
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Figure 5. Effect of doping on the critical field h c . In the left panel we show the effect 
of S on the SF-P transition, close to the A = case. The right panel shows the value 
of the critical field h c for the AF-SF transition for two values of A. The t, t', and J 
parameters are those of Table ^ 

finite range of A values). On the other hand, the effect for values such that A ~ 1 is to 
introduce an increase of h c . Although our calculation cannot be extended quantitatively 
to La2- x Sr x Cu04 mainly because the spin flop transition in this material is due to the 
Dzyaloshinskii-Moriya interaction we do not expect qualitative changes in what respects 
the behavior of h c with doping. [26j That is, it is possible to account at the mean field 
level for the decrease of h c with doping, independently of the details of the interaction, 
as lone AF-SF transition exists. 

Let us see now how the zero-temperature picture evolves when we extend the 
analysis to finite temperatures. Since the Bose-Einstein condensation temperature is 
of the order of the doping, Tbe ~ St, we can draw a finite temperature phase diagram 
using the same formulation we used for zero T as long as T < T BE . We have stayed 
in this regime. As for the pure Heisenberg case, we obtain a bicritical point where the 
first order AF-SF transition line meet the two second order lines, describing AF-P and 
the SF-P transitions. Comparing with the pure case, the main changes are: (i) the 
bicritical point (T&, hb) moves to lower temperatures; (ii) the SF region shrinks as we 
dope the system up to a point where only the AF and the P phases remain (the area of 
the AF zone is also reduced). The disappearance of the SF phase with doping before the 
AF phase, is related to the zero temperature dependence of the average value S upon 
doping that is seen in figure 01 where the S attains its paramagnetic value before the Sa 
and Sb do so. The aspects discussed above are illuminated in figure El where in the left 
panel the phase diagram is depicted, and in the right one we plot the dependence of the 
average values of S, Sa, and Sb, and the angle cf> in the SF phase along the first order 
AF-SF transition line. We see that close to the bicritical point, and in a reduced range 
of temperatures, the angle (p has a fast variation. Therefore, for an experiment probing 
the AF-SF transition with temperature, at different fields close to the bicritical point 
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Figure 6. Finite temperature phase diagram of the 3D t — t' — J model, for A=2 and 
1 — S = n = 0.92. The t, t', and J parameters are those of Table The left panel 
shows the phase diagram, and the right panels show S, Sa, Sb, and 4> as function of 
T along the AF to SP first-order transition line. 



we may see a large or small jump in the magnetization, depending on the field strength. 
Such a behavior could be easily observed performing magnetization measurements of 
the type presented for L^-^Sr^CuCU, jTH] but with the magnetic field applied along 
the Cu04 planes. Again, our discussion of the relation between our findings and the 
experimental results in La 2 _ :E Sr :r Cu04 applies with the limitations discussed above. 

4. Characterization of the spin flop phase 

Let us characterize the metallic SF phase. We consider both magnetic and charge 
transport properties. We first compute the magnetic longitudinal and transverse 
susceptibility. After we study the optical conductivity, computing both the Drude weight 
and the regular part of the conductivity. 

4-1. Magnetic susceptibility 

The dynamic magnetic susceptibility, Xa,/3(Q,iuj n ), is defined as 

XaA^n) = f dre^iTS^q^Spi-q^)) , (13) 
J o 

where 

%) = ^E fl +q ,atf' P) fk,«- (14) 
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In the SF phase we can define four types of Green's functions, generically written 

as 

G a Ak,p,r) = -(Tf k , a (r)fU°)) ■ ( 15 ) 
For each case the results in the Matsubara representation are 

8k, P [iv n -ei(k- Q)} 



G TiT (A;,p, iu n ) 
G[,\{k,p, iuj n ) 



[iuj n -e±(k- Q)][iuJ n - t\{k)\ - A 2 
AG TiT (fc + Q,p,iu n ) 
iun-e^k) 

h, P l iuJ n - e T (fc + Q)] 

[iu n - ei(k)][iu n - e T (A; + Q)] - A 2 
AG U (& - Q,P,iv n ) 



iu n -e^(k) 

Since G^y(k,p,iu n ) and G^i(k,p,iuj n ) are diagonal in momentum space, the only 
non-zero off-diagonal Green's functions in spin space are G^{p — Q,p,iu n ) and 
G^i(p + Q,p,iuJ n )- 

At the mean field level, the longitudinal magnetic susceptibility, Xz,z.i an d the 
transverse magnetic susceptibility, X-,+> are obtained using the definition (|T5|l combined 
with the Green's functions (115)1 . The general form of the susceptibility is 

1 2 

iV P i,l=l 

x[M jM ^(p,q)] 2 , (16) 



where 



and 



T[E j (p + q),E l (p),iu n ] = \ r-—, 17 

iuj n + Ejip + q) - Ei(p) 

M h i, z M q) = \[RiAp + q)RM - r iAp + q)RiM • ( 18 ) 

Mj,i,-,+(p, q) = R Uj (p + q)Ru(P ~ Q) > (19) 



Here f(x) = (1 + e x / T ) 1 and the R a j factors are given in |Appendix A 



The behavior of Xz,z(®-> 0) anc ^ X-,+ (0,0) is presented in figure (J7J). Above the 
SF-P transition line defined by the phase diagram of figure El the behavior of the 
susceptibility is that of a paramagnet in a magnetic field. As h — > 0, x+- — %Xzz- 
Below the transition temperature, there is a sudden drop in the magnetic susceptibility, 
followed by the same abrupt behavior of the SP angle <fi. As the temperature is further 
reduced the susceptibility shows an up turn at low temperature, which is typical of the 
susceptibility of metallic antiferromagnets, as is the case of PtsFe. j2Z| 
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Figure 7. Spin susceptibility Xz,z(0,0) and x-,+ (0,0) at finite temperature in the 
spin-flop region of the phase diagram for A=2, h = 0.6 and n = 0.92. The t, t' , and J 
parameters are those of Table ^ 



4-2. Optical conductivity 

The response of the system to an electromagnetic field is obtained from the optical 
conductivity. This quantity in defined as 



\ 1 (K xx ) +A xx (q,u) 

where A xx (q,u) is the retarded current-current correlation function, obtained from the 
corresponding Matsubara correlation function 

A xx (q,iu n ) = [ P dre^(T T f x (q,r)f x (-q,0)) , (21) 



and K xx is given by 

K xx = £e(fc + Q)f\ MQ hk +Q + e(k)fl k f l>k . (22) 

k 

Here e{k) = —2t8 cosq x — 4t'S cosg^cosg,, + cosg z ) and j x (q) is a Fourier component of 
the current operator. [2E] Writing cr(q, uj) = a'(q,u) + ia" (q, u) , that is, separating the 
real and imaginary parts, the real part reads 

AM = - ^s(u)[(k xx ) + KM")\ + 

= n5(uj)D + cr reg (q,uj) , (23) 
where D is the charge stiffness or Drude weight, given by 

D = ~[(K„)+h!JPM- (24) 
The zero momentum conductivity is given by 

A£,(0,u;) 



a rP „(uj 



nJ2 E M md S(u + Ej-E m ), (25) 
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Figure 8. Optical conductivity and Drude weight of the 3D t — t' — J model, for A=2 
and n = 0.92, h=0.6. The t, t', and J parameters are those of Tabled I n the left 
panels we show the effect of the doping on the zero temperature Drude weight D and 
the effect of the temperature on D as the SF-P transition line is crossed. In the right 
panel, the regular part a reg of the optical conductivity is depicted. 

with 

M m>j = \f(E 3 ) - f{E m )]\j 2 (p + Q)Rl m Rl j 
+ 2j(p + Q)j{p) + R^ m R 1d R lm R u 

+ f{p)Rl m Rl 3 \ , (26) 

where j(k) = 2t8 sin q x + At' 8 sin q x (cos q y + cos q z ). 

A finite Drude weight establishes an infinite d.c. conductivity. If the system is an 
insulator D is zero. On the other hand, a reg (uj), establishes the absorption of finite 
frequency light by the system. In figure |H1 we show the effect of the doping and of the 
temperature on D, as well as the effect of the doping on a re g{oj) at zero temperature. 
At zero temperature we realize that the effect of increasing 5 is two fold: (i) it increases 
D (more carriers available); (ii) it shifts a reg {uj) to lower values in frequency, because 
the gap decreases with 5. When 5^0, both D and a reg {uj) vanish. We expect this 
result to hold beyond our mean field analysis. 29J At finite temperature and constant 5, 
there are two situations. Above the SF-P transition the Drude weight exhausts the sum 
rule obeyed by c(uj), and therefore there is no finite energy absorption. When the SF-P 
transition line is crossed, there is a reduction of D, with a transfer of spectral weight to 
finite energies. The diminishing of D is quite abrupt in a small range of temperatures, 
and this fact opens the possibility of having a large magneto-resistance in these systems, 
when scattering from impurities is included. 
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5. Conclusions 

In conclusion, we have, in this article, described the mean field theory of spin-flop 
transitions in a doped antiferromagnets. Further we have characterized the spin-flop 
phase from the point of view of its magnetic and charge-transport properties. We have 
shown the AF phase is characterized by a field-induced ferrimagnetic ground state, at 
odds with the pure Heisenberg case, and that a SF phase exists in a finite range of 
A and h. The Drude weight of the system was shown to present an abrupt decrease 
in a reduced range of temperatures when the system enters the spin-flop phase, and 
a similar behavior occurs with the magnetic susceptibility. We have also studied the 
behavior of the critical lines with the doping, and found that depending on the value of 
A the critical field associated with the AF to SF transition-line may increase (A ~ 1) or 
decrease (A > 1) with the doping 5. We have also shown that the behavior of the critical 
field associated with the SF to P transition is non monotonous with 5. Although, the 
full magnetic behavior of doped La2Cu04 in a magnetic field cannot be quantitatively 
described by an antiferromagnet with an anisotropic exchange, since the origin of the 
spin-flop phase in this material is due to the Dzyaloshinskii-Moriya interaction, we 
believe that some of the above qualitative findings still apply. The main reason, we 
believe, is the fact that some physical properties will not depend so much on the details 
of the interaction but on the existence of a physical mechanism allowing for a spin flop 
transition. Only this reasoning may explain why the upper right panel of Fig. 6 exhibits 
a decreasing of the lower critical field with doping, as observed in doped La 2 Cu0 4 . Also 
the existence of an incommensurate spin flop phase will not depend on the detail of the 
interaction but only on the competition between the kinetic and interaction energies. 
The reason why an incommensurate spin flop state is not observed in pure La2Cu04 is 
precisely the absence of the kinetic energy term. On the contrary, details as the line 
borders separating the different phases, the values of the incommensurate momentum 
and the amount of reduction of the critical field with doping will certainly be interaction 
dependent. 
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Appendix A. Green's function decomposition, mean field equations and 
susceptibility 

In the spin flop phase all the Green's functions can be cast in the form 

G a ,,(p,^„) = £^%. (A.l) 
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In particular we have that 

G u (p + Q,iu n ) =E [R ' AP)? 



GuiP ^ n) ~ g iu> n - E s {p) ' 



where the coherence factors are 

[^t,i(p)] 2 = [^2(p)] 2 = -(:i + 



2 

E[^>)] 2 , 



i2 A> i(p)^T,i(p) = -^,2(p)i2T,a(p) = ~5= » ( A -2) 

with a(p) = [A(p)] 2 + 4A 2 and A(p) = e T (p + Q) - e^p). 

The mean field equations can be expressed in terms of the coherence factors as 

iV fc 3=1 

Ssin(20)[Az- 7 (g)] = -7(Q)cos^X)E 2i2 T^w/(^) + 

A *w-^EE(4. - #iM E >) » ( A - 4 ) 

fc jr' = l 

%A cos 2 + 7 (g) sin 2 0] = Az cos ]T £(7^. - i? 2 ,,)/^) + 



2Ar * ;=i 



1 2 



7 (g) sin 0— EE IR^RuttEj) • (A.5) 

The above equations hold for both the commensurate and incommensurate cases. In 
the incommensurate case an additional equation is required for the determination of Qj 
(j = x, y, z)- the incommensurate momentum value. This equation reads 

+ ^ ^ * sin * = • < A - 6 > 

the partial derivatives of the quasi-particle dispersions E a are straightforward but give 
long equations we omit here. 
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The full expression for the Xzz(0, 0) susceptibility useful for any value of the doping 
5 is given by 

2A 2 

x»(o,o)= -JC^/^j^p 

1 [A(fc)] 2 /(E a )/(-E a ) 
+ ,i ± 4 Ta(*) ' (A - 7) 
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